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We present improved measurements of CP- violation parameters in the decays B° — > ir + ir~ , B° — > 
K + n-, and B° -> ttV, and of the branching fractions for B° -> tv°-k and B° -> K°tv°. The results 
are obtained with the full data set collected at the T(4S) resonance by the BABAR experiment at the 
PEP-II asymmetric-energy B factory at the SLAC National Accelerator Laboratory, corresponding 
to 467 ± 5 million BB pairs. We find the CP-violation parameter values and branching fractions 

S w+n - = -0.68 ± 0.10 ± 0.03, 
C K+V - = -0.25 ± 0.08 ± 0.02, 
A K -„+ = -O.107±O-016±g;S8S, 
C n o„o = -0.43 ± 0.26 ± 0.05, 
B(B° -> ttV ) = (1.83 ± 0.21 ± 0.13) x 10 -6 , 
B(B° -> A'V ) = (10.1 ± 0.6 ± 0.4) x 10 



(i 



where in each case, the first uncertainties are statistical and the second are systematic. We observe 
CP violation with a significance of 6.7 standard deviations for B° — > -k + -k~ and 6.1 standard devi- 
ations for B° — > K + iv~ , including systematic uncertainties. Constraints on the Unitarity Triangle 
angle a are determined from the isospin relations among the B — > tt-k rates and asymmetries. Con- 
sidering only the solution preferred by the Standard Model, we find a to be in the range [71°, 109°] 
at the 68% confidence level. 
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I. INTRODUCTION 
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Large CP-violating effects [lj in the P-meson sys- 
tem are among the most remarkable predictions of 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing 
model Q. These predictions have been confirmed by 
the BABAR and Belle Collaborations, most precisely in 
b — s- ccs decays of P° mesons to CP eigenstates [H, |j] . 

Effective constraints on physics beyond the Standard 
Model (SM) are provided by high-precision measure- 
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ments of quantities whose SM predictions are subject to 
only small theoretical uncertainties. Many experimental 
and theoretical uncertainties partially cancel in the cal- 
culation of CP- violating asymmetries. This makes CP- 
violation measurements a sensitive probe for effects of 
yet-undiscovered additional interactions and heavy par- 
ticles that are introduced by extensions to the SM. All 
measurements of CP violation to date, including those 
involving the decay modes studied here , are m 
agreement with the indirect predictions from global SM 
fits [l(J El, which are based on measurements of the 
magnitudes of the elements Vij of the CKM quark-mixing 
matrix. This strongly constrains 1 1 21 the flavor structure 
of SM extensions. 

The CKM-matrix unitarity-triangle angle a = 
arg [— ^td^tb/Kid^ub] i s measured through interference 
between two decay amplitudes, where one amplitude in- 
volves B°-B° mixing. Multiple measurements of a, with 
different decays, further test the consistency of the CKM 
model. The time-dependent asymmetry in B° — > ir + ir~ 
decays is proportional to sin2a in the limit that only the 
b — > u ( "tree" ) quark-level amplitude contributes to this 
decay. In the presence of b — > d ("penguin") amplitudes, 
the time-dependent asymmetry in B° — > n + ir~ is modi- 
fied to 



a(At) 



\A(At)\ 2 \A(At)f 



\A(At)\* + \A(At)\2 



sin (Am^At) — C\-i 



cos (ArndAt), 
(1) 



where At is the difference between the proper decay times 



of the B meson that undergoes the B 



decay 



(the signal B) and the other B meson in the event (the 
tag B), Arrid is the B°-B° mixing frequency, A is the 
B° —> 7T + 7r~ decay amplitude, A is the CP-conjugate 
amplitude, and 



C- 



\A\ 2 -\Af 
|A| 2 + L4| 2 ' 



l-C2 +w _sin(2a-2Aa^). (2) 



Both the direct CP asymmetry C v + V - and the phase 
Ao^,,. may differ from zero due to the penguin contribu- 
tion to the decay amplitudes. 

The magnitude and relative phase of the penguin con- 
tribution to the asymmetry S n + n - may be determined 
with an analysis of isospin relations between the B — > tttt 
decay amplitudes [llj]. The amplitudes A^ of the B — s- 
ir 1 ^ decays and A y of the B — > 7TV- 7 decays satisfy the 
relations 

A+° = -La+-+A qo , 
\/2 



A~ u = —^=A 

V2 



A 



mi 



(3) 



The shapes of the triangles corresponding to these 
isospin relations are determined from measurements of 



the branching fractions and time-integrated CP asym- 
metries for each of the B — > tttt decays. Gluonic penguin 
amplitudes do not contribute to the AI = 3/2 decay 
£? ± — > 7r^~ 7r® . Therefore, neglecting electroweak (EW) 
penguin amplitudes, the amplitudes A +0 and A~° are 
equal. From the different shapes of the triangles for the 
B and B decay amplitudes, a constraint on Aa„ can be 
determined to within a four-fold ambiguity. 

The phenomenology of the B — > irir system has been 
thoroughly studied in a number of theoretical frame- 
works and models [l4j]. Predictions for the relative size 
and phase of the penguin contribution vary considerably. 
Therefore, increasingly precise measurements will help 
distinguish among different theoretical approaches and 
add to our understanding of hadronic B decays. 

The measured rates and direct CP-violating asymme- 
tries in B ->■ Kit decays [f| 0, H, ITsT - fl^ l reveal puz- 
zling features that could indicate significant contribu- 
tions from EW penguin amplitudes [13, |20j. Various 
methods have been proposed for isolating the SM con- 
tribution to this process in order to test for signs of new 
physics. This includes sum rules derived from /7-spin 
symmetry, which relate the rates and asymmetries for 
the decays of charged or neutral B mesons to K + u~ , 
K+n°, K°tt°, and K n+ [IHH, and SU(3) symmetry, 
used to make predictions for the Kn system based on 
hadronic parameters extracted from the tttt system [191 ]. 

This article is organized as follows. The BABAR de- 
tector and the data used in these measurements are de- 
scribed in Section[TTl In Scction lTOl we outline the analysis 
method, including the event selection and the fits used 
to extract the parameters of interest. The results of the 
data analysis are given in Section llVl The extraction of a 
and Aa^Tr is described in Section [Vj and we summarize 
in Section IVT1 



II. THE BABAR DETECTOR AND DATA SET 

In the BABAR detector [23|, charged particles are de- 
tected and their momenta are measured by the combi- 
nation of a five-layer double-sided silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) that cov- 
ers 92% of the solid angle in the T(4S) center-of-mass 
(cm.) frame, both operating in a 1.5 T uniform magnetic 
field. Discrimination between charged pions, kaons, and 
protons is obtained from ionization (dE/dx) measure- 
ments in the DCH and from an internally reflecting ring- 
imaging Cherenkov detector (DIRC), which covers 84% 
of the cm. solid angle in the central region of the BABAR 
detector and has a 91% reconstruction efficiency for pions 
and kaons with momenta above 1.5 GeV/c. Photons and 
electrons are identified and their energies are measured 
with an electromagnetic calorimeter (EMC) consisting 
of 6580 CsI(Tl) crystals. The photon energy resolution 
is a E /E = |2.3/£:(GeV) 1 / 4 ® 1.4}%, and the photon 
angular resolution relative to the interaction point is 
a e = 4.16/v/P(GeV) mrad HJ. 
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The data used in this analysis were collected during 
the period 1999-2007 with the BABAR detector at the 
PEP-II asymmetric-energy B-meson factory at the SLAC 
National Accelerator Laboratory. A total of 467 ± 5 mil- 
lion BB pairs were used. Relative to previous EABAR 
measurements 043) roughly 22% more BB pairs have 
been added to the analyzed data set, and improvements 
have been introduced to the analysis technique, boosting 
the signal significance. These improvements include bet- 
ter reconstruction of charged-particle tracks, improved 
hadron-identification and flavor-tagging algorithms, and 
optimal selection of tracks and calorimeter clusters for 
calculation of event-shape variables. 

Samples of Monte Carlo (MC) simulated events are an- 
alyzed with the same reconstruction and analysis proce- 
dures as used for the data, following a GEANT4-based (2f| 
detailed detector simulation (23| . The MC samples in- 
clude e + e~ — > qq continuum background events gener- 
ated with JETSET [51 and T(4S) ->• BB decays gen- 
erated with EvtGen [27| and JETSET, including both 
signal and background B-meson decays. 



III. EVENT SELECTION AND ANALYSIS 
METHOD 

Many elements of the measurements discussed in this 
paper are common to the decay modes [28[ B° — > h + h'~ 
(where h,h' = tt or K), B° -> 7r°7r°, and B° K°tt°. 
The signal i?-meson candidates (B ICC ) are formed by 
combining two particles, each of which is a charged- 
particle track, a ir° candidate, or a K° candidate. The 
event selection differs for each mode, and is described 
below. 

The number of B decays and the corresponding CP 
asymmetries are determined with extended unbinned 
maximum likelihood (ML) fits to variables described be- 
low. The likelihood is given by the expression 



M 



N 



C = exp 



n 



(4) 



where N is the number of events, the sums are over the 
event categories M, rii is the event yield for each category 
as described below, and the probability-density function 
(PDF) Vi describes the distribution of the variables Xj in 
terms of parameters on. The PDF functional forms are 
discussed in Sections IIII CI and IIII Dl 



A. Track and K° Selection 

In the B° — » h + h'~ mode, we require charged-particle 
tracks to have at least 12 DCH hits and to lie in the 
polar-angle region 0.35 < 9 < 2.40 with respect to the 
beam direction. The track impact parameter relative to 
the e + e~ collision axis must be smaller than 1.5 cm in 



the plane perpendicular to the beam axis and 2.5 cm in 
the direction along the axis. 

In order for DIRC information to be used for particle 
identification, we require that each track have its asso- 
ciated Cherenkov angle (9c) measured with at least six 
Cherenkov photons, where the value of 9c is required 
to be within 4.0 standard deviations (a) of either the 
pion or kaon hypothesis. This removes candidates con- 
taining a high-momentum proton. Tracks from electrons 
are removed based primarily on a comparison of the 
track momentum and the associated energy deposition in 
the EMC, with additional information provided by DCH 
dE/dx and DIRC 9c measurements. 

The ionization energy loss in the DCH is used either in 
combination with DIRC information or alone. This leads 
to a 35% increase in the B° — > h + h!~ reconstruction effi- 
ciency relative to the use of only tracks with good DIRC 
information. A detailed DCH dE/dx calibration devel- 
oped for the B° — > h + h!~ analysis takes into account 
variations in the mean and resolution of dE/dx measure- 
ment values with respect to changes in the DCH running 
conditions over time, as well as the track's charge, po- 
lar and azimuthal angles, and number of ionization sam- 
ples. The calibration is performed with large high-purity 
samples (> 10 6 events) of protons from A — > p7r~, pi- 
ons and kaons from D* + — > D a ir + (D° — > K~tt + ), and 
K® — > 7r + 7r~ decays that occur in the vicinity of the in- 
teraction region. 



Candidates for the decay Kg 



7r n 7r are recon- 



structed from pairs of oppositely-charged tracks. The 
two-track combinations are required to form a vertex 
with a x 2 probability greater than 0.001 and a 7r + 7r~ in- 
variant mass within 11.2 MeV/c 2 , corresponding to 3.7<r, 



of the nominal K$ mass |2£ 



B. 7T° Selection 

We form tt° — > 77 candidates from pairs of clusters in 
the EMC that are isolated from any charged track. Clus- 
ters are required to have a lateral profile of energy de- 
position consistent with that of a photon and to have an 
energy £ 7 > 30 MeV for B° 7r°7r° and E 1 > 50 MeV 
for B° — > K®tt°. We require ir° candidates to lie in the 
invariant-mass range 110 < to 77 < 160 MeV/c 2 . 

For the B° — > 7r°7r° mode, we also use 7r° candidates 
from a single EMC cluster containing two adjacent pho- 
tons (a merged 7r°), or one EMC cluster and two tracks 
from a photon conversion to an e + e~ pair inside the de- 
tector. To reduce the background from random photon 
combinations, the angle 9 y between the photon momen- 
tum vector in the tt° rest frame and the 7r° momen- 
tum vector in the laboratory frame is required to sat- 
isfy |cos6* 7 | < 0.95. The ir° candidates are fitted kine- 
matically with their mass constrained to the nominal tt° 
mass [29l | . 

Photon conversions are selected from pairs of 
oppositely-charged electron-candidate tracks with an in- 
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variant mass below 30 MeV/c 2 whose combined momen- 
tum vector points away from the beam spot. The con- 
version point is required to lie within detector material 
layers. Converted photons are combined with photons 
from single EMC clusters to form ir° candidates. 

Single EMC clusters containing two photons are se- 
lected with the transverse second moment, S = J^. Ei x 
(Acti) 2 /E, where Ei is the energy in each CsI(Tl) crys- 
tal and Acti is the angle between the cluster centroid 
and the crystal. The second moment is used to distin- 
guish merged ir° candidates from both single photons and 
neutral hadrons. 



C. B° -> TT+TT", B° ->■ K+n~, and B° -> 7r°7r° 

Two kinematic variables are used in the B° —> h + h'~ 
and B° — > tt q tt° analyses to separate B- meson decays 
from the large e + e~ — > qq (q — u, d, s, c) combina- 
toric background (Hf. One variable is the beam-energy- 
substituted mass toes = \J (s/2 + • pb) 2 /E 2 — p 2 B , 
where y 7 ? is the total e + e~ cm. energy, (£^,pj) is the 
four-momentum of the initial e + e~ system in the labora- 
tory frame, and ps is the laboratory momentum of the 
B candidate. The second variable is AE = E* B — y/s/2, 
where E* B is the energy of the B candidate in the cm. 
frame. 

To further separate B decays from the qq background, 
we use two additional topological variables that take 
advantage of the two-jet nature of qq events and the 
isotropic particle distribution of e + e~ — > BB events. 
The first variable is the absolute value of the cosine of 
the angle 9 S between the sphericity axis [30[ of the decay 
products of the B candidate and the sphericity axis of the 
remaining tracks and neutral clusters in the event, com- 
puted in the cm. frame. The distribution of this variable 
peaks at 1 for the jet-like qq events and is uniform for B 
decays. We require | cos# s | < 0.91 for B° — »• h + h!~ and 
Icos^sl < 0.7 for B° —> 7r°7r°, where a tighter require- 
ment is needed due to the higher background. For the 
B° — > h + h'~ mode, we remove a small remaining back- 
ground from e + e~ — > r + r~ events by further requi ring 
that the normalized second Fox- Wolfram moment [3l| 
satisfy R 2 < 0.7. 

To improve the discrimination against qq events, a 
Fisher discriminant J- is formed as a linear combination 
of the sums Lq = £ilP*l and L\ = |p*| cos 2 0?, 
where p* are the momenta and 9* are the angles with re- 
spect to the thrust axis [12] of the B candidate, both in 
the cm. frame, of all tracks and clusters not used to re- 
construct the signal £?-meson candidate. The T variable 
takes advantage of the fact that much of the momentum 
flow in qq events is along the thrust axis. In the case of 
B° — » 7r°7r°, we improve the sensitivity to signal events 
by combining T with three other event-shape variables 
in a neural network. The first variables is | cos 6 S | , de- 
scribed above. The second is |cos0g|, where 9* B is the 
angle between the momentum vector of the signal B and 



the beam axis. The | cos6>g| distribution of qq events is 
uniform, while that of signal events is proportional to 
sin 2 9* B . The third variable is | cos 9* T | , where 9* T is the 
angle between the thrust axis of the signal i?-meson's 
daughters and the beam axis. Both 9* B and 9^ are cal- 
culated in the cm. frame. The characteristics of the 
| cos 9?p | distributions are similar to those of | cos 9 S | . 

1. B° ->■ TT+TT" and B° ->• K + tt~ 

We reconstruct the candidate decays B rcc — > h + h'~ 
from pairs of oppositely-charged tracks that are consis- 
tent with originating from a common decay point with 
a x 2 probability of at least 0.001. The remaining par- 
ticles are examined to infer whether the other B meson 
in the event (-Stag) decayed as a B° or B° (flavor tag). 
We perform an unbinned extended ML fit to separate 
B° — > 7r + 7r~ and B° — > K + ix~ decays and determine 
simultaneously their CP-violating asymmetries S w + W -, 
Cw+yj.— , and 

_ B(B -> K-TT+) - B(B -> K+TT-) 
Ak ^ + ~ BjB K-TT+) + B(B ->■ K+TT-) ' ^ ' 

as well as the signal and background yields and PDF 
parameters. The fit uses 9c, dE/dx, AE, m ES , J 7 , B tag 
flavor, and At information. 

The value of AE is calculated assuming that both 
tracks are charged pions. The B° — s- n + ir~ signal is 
described by a Gaussian distribution for AE, with a res- 
olution of 29 MeV. For each kaon in the final state, the 
AE peak position is shifted from zero by an amount 
that depends on the kaon momentum, with an average 
shift of -45 MeV. We require \AE\ < 0.150 GeV. The 
wide range in AE allows us to separate B° decays to 
the four final states ir + ir~ , K + ir~, ir + K~, and K + K~ 
in a single fit. The analysis is not optimized for mea- 
suring the K + K~ final state, which is treated as back- 
ground. The toes resolution is 2.6 MeV/c 2 . We require 
toes > 5.20 GeV/c 2 , with events in the large range below 
the signal peak allowing the fit to effectively determine 
the background shape parameters. 

We construct 9c PDFs for the pion and kaon hypothe- 
ses, and dE/dx PDFs for the pion, kaon, and proton 
hypotheses, separately for each charge. The K-n sepa- 
rations provided by 9c and dE/dx are complementary: 
for 9c, the separation varies from 2.5cr at 4.5 GeV/c to 
13a at 1.5 GeV/c, while for dE/dx it varies from less than 
1.0a at 1.5 GeV/c to 1.9a at 4.5 GeV/c (Fig.HJ. For more 
details, see Ref. [5[. 

We use a multivariate technique [33| to determine the 
flavor of the -Btag- Separate neural networks are trained 
to identify leptons from B decays, kaons from D decays, 
and soft pions from D* decays. Events are assigned to 
one of seven mutually exclusive tagging categories (one 
category being untagged events) based on the estimated 
average mistag probability and the source of the tag- 
ging information. The quality of tagging is expressed 
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in terms of the effective efficiency Q = J2k e k(^ ~ % w k) 2 , 
where and w^ are the efficiencies and mistag prob- 
abilities, respectively, for events tagged in category k. 
The difference between the mistag probabilities for B° 
and B° mesons is given by Aw = wgo — wg . Table U 
summarizes the tagging performance measured in a large 
data sample of fully-reconstructed neutral B^ av decays 



to £>M- (tt+ 



a+) $4 



FIG. 1: The average expected K—n separation, in units of un- 
certainty, provided by the DIRC angle 6c and DCH dE/dx for 
kaons and pions from B° —¥ K + ty~ decays in the laboratory- 
frame polar angle range 0.35 < 6 < 2.40, as a function of 
laboratory-frame momentum. 



TABLE I: Average tagging efficiency e, average mistag frac- 
tion w, mistag fraction difference Aw = w(B°) — w(B°), and 
effective tagging efficiency Q for signal events in each tagging 
category (except the untagged category). 



Category 



:(%) 



w (%) Aw (%) 



•(%) 



Lepton 
Kaon I 
Kaon II 
KaonPion 13.67 
Pion 
Other 



Total 



8.96 ± 0.07 2.9 ± 0.3 

10.81 ± 0.07 5.3 ± 0.3 

17.18 ± 0.09 14.5 ± 0.3 
0.08 23.3 ± 0.4 

14.19 ± 0.08 32.6 ± 0.4 
9.55 ± 0.07 41.5 ± 0.5 



0.2 ± 0.5 7.95 ± 0.11 
0.0 ± 0.6 8.64 ± 0.14 
0.4 ± 0.6 8.64 ± 0.17 
-0.6 ± 0.7 3.91 ± 0.12 
5.1 ± 0.7 1.73 ± 0.09 
3.8 ± 0.8 0.28 ± 0.04 



31.1 ± 0.3 



The time difference At — Az//3"/c is obtained from 
the known boost of the e + e~ system (/?7 = 0.56) and 
the measured distance Az along the beam (z) axis be- 
tween the B rec and -Btag decay vertices. A description of 
the inclusive reconstruction of the Stag vertex is given 
in Ref. [35j. We require |At| < 20 ps and (JAt < 2.5 ps, 
where cta* is the uncertainty on At, estimated separately 
for each event. The signal At PDF for B° —> tt + tt~ is 
given by 



ff~ (At mcas ) 



,-|At|/ T 

At 



{(IT Aw) ±(l-2u>fc) sin (ArridAt) — C^+^- cos (Am d At)] } ® R(At 

mcas 

At), (6) 



r 



where f£ (fj7) indicates a B° (B°) flavor tag and the 
index k indicates the tagging category. The resolution 
function i?(At moas — At) for signal candidates is a sum of 
three Gaussian functions, identical to the one described 
in Ref. |35| . with parameters determined from a fit to 
the -Bfl a v sample, which includes events in all seven tag- 
ging categories. The background At distribution is mod- 
eled as the sum of three Gaussians, with parameters, 
common for all tagging categories, determined simulta- 
neously with the CP violation parameters in the ML fit 
to the B lec — > h + h!~ sample. 

The ML fit PDF includes 28 components. Of these, 
24 components correspond to B° signal decays and back- 
ground events with the final states 7r + 7r~, K + ir~, K~ir + , 
and K + K~, where either the positively-charged track, 
the negatively-charged track, or both have good DIRC 



information (2 x 4 x 3 — 24 components). Four addi- 
tional components correspond to jot - , pK~, n + p and 
K + p background events, where the (anti)proton has no 
DIRC information. The K ir^ event yields riK±-nT are 
parameterized in terms of the asymmetry A T 3^ + and av- 
erage yield n K7r as n K±7rT = n Klx (l ^fA T ^l w+ ) /2. All 
other event yields are products of the fraction of events in 
each tagging category, taken from i?fl av events, and the 
total event yield. The background PDFs are a thresh- 
old function (3(| for toes and a second-order polyno- 
mial for AE. The T PDF is a sum of two asymmetric 
Gaussians for both signal and background. We use large 
samples of simulated B decays to investigate the effects 
of backgrounds from other B decays on the determina- 
tion of the CP-violating asymmetries in B° — > ir + ir~ and 
B° —> K + ir~ , and find them to be negligible. 
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2. B° -> 7r°7r° 

B° — > 7T 7r° events are identified with an ML fit to the 
variables Toes, AE, and the output AW of the event- 
shape neural network. We require toes > 5-20 GeV/c 2 
and \AE\ < 0.2 GeV. Since tails in the EMC response 
produce a correlation between Toes an d AE, a two- 
dimensional binned PDF, derived from the signal MC 
sample, is used to describe signal PDF. The AW distribu- 
tion is divided into ten bins (with each bin approximately 
equally populated by signal events) and described by a 
nine-bin step-function PDF with values taken from the 
MC and fixed in the fit. -Bfl av data are used to verify 
that the MC accurately reproduces the AW distribution. 
The qq background PDFs are a threshold function [36[ 
for TOeSj a second-order polynomial for AE, and a para- 
metric step function for AW. For qq events, AW is not 
distributed uniformly across the bins but rises sharply 
toward the highest bins. We see a small correlation of 
2.5% between the shape parameter of the Toes threshold 
function and the AW bin number, and this relation is 
taken into account in the fit. All qq background PDF- 
parameter values are determined by the ML fit. 

The decays B+ ->■ p+Tr and B° K°-k° (K° 
7r°7r°) add 71 ± 10 background events to B° -t 7r°7r° 
and are included as an additional component in the ML 
fit. We model these B-decay background events with a 
two-dimensional binned PDF in toes and AE, and with 
a step function for AW. The shapes of these PDFs are 
taken from MC simulation, and their event yields and 
asymmetries are fixed in the fit and are later varied to 
evaluate systematic uncertainties. 

The time-integrated CP asymmetry is measured by the 
-B-flavor tagging algorithm described above. The fraction 
of events in each tagging category is constrained to the 
corresponding fraction determined from MC simulation. 
The PDF event yields for the B° — > tt°tt° signal are given 
by the expression 

n^o )k = -fkN^o l-s-,(l-2x)(l-2w fe )C 7r o 7r o , (7) 

where is the fraction of events in tagging category k, 
N n o w o is the number of B° — > TT n candidate decays, \ 
is the time-integrated B° mixing probability [29], Sj = 
+1(-1) when the B tag is a B° (B°), and 



|^00|2_ |^00|! 



(8) 



|^00|2 + |^00|2 

is the direct CP asymmetry in B° -> 7r°7r°. 

D. B" -s. AT°7r 

CP-violation parameters for B° — > K®tt° have been 
reported in Ref. 0. Here we describe the measurement 
of the branching fraction for this mode. 



For each B° — > K®ir° candidate, two independent kine- 
matic variables are computed. The first variable is the 
invariant mass me of the B vcc . The second variable is 
the invariant (missing) mass TO m ; ss of the -Btag , computed 
from the magnitude of the difference between the four- 
momentum of the initial e + e" system and that of the 
B roc , after applying a -B°-mass constraint to the B lec [37| . 
For signal decays, raj and TO m i ss peak near the B° mass 
with resolutions of about 36 and 5.3 MeV/c 2 , respectively. 
Since the linear correlation coefficient between mg and 
TO m i SS vanishes, these variables yield better separation of 
signal from background than toes and AE. Both the 
tub and m m iss distributions exhibit a low-side tail due 
to leakage of energy out of the EMC. We select candi- 
dates within the ranges 5.13 < tub < 5.43 GeV/c 2 and 
5.11 < TO m i ss < 5.31 GeV/c 2 , which include a signal peak 
and a "sideband" region for background characterization. 
In events with more than one reconstructed candidate 
(0.8% of the total), we select the candidate with the 
smallest \ 2 — X«=7r°,if ( TO i ~ m '%)' 2 / m i > where m l (to-) 
is the measured (nominal) mass and a mi is the estimated 
uncertainty on the measured mass of particle i. 

We exploit topological observables, computed in the 
cm. frame, to discriminate jet-like e + e~ — > qq events 
from the nearly spherical BB events. In order to reduce 
the number of background events, we require -L2/L0 < 
0.55, where Lj = Ip^I cos j 6** and 0* are computed 
with respect to the sphericity axis [3(| of the B ICC can- 
didate. Taking advantage of the fact that signal events 
follow a 1 — cos 2 6* B distribution while the background 
is flat, we select events with | cos 9* B \ < 0.9. Using a 
full detector simulation, we estimate that our selection 
retains (34.2 ± 1.2)% of the signal events, where the un- 
certainty includes both statistical and systematic con- 
tributions. The selected sample of B° —> K®ir° candi- 
dates is dominated by random K®tt° combinations from 
e + e~ — > qq events. Using large samples of simulated BB 
events, we find that backgrounds from other -B-meson 
decays are small, of order 0.1%. Therefore, this type of 
background is not included in the fit described below, 
and this is accounted for in the evaluation of systematic 
uncertainties (see Section HV C| . 

We extract the signal yield from an extended unbinned 
ML fit to to b j TO m iss, L2/L0, cos9 Bl the flavor-tag, and 
the decay time and its error. The use of tagging and 
decay-time information in the ML fit further improves 
discrimination between signal and background. Since in 
the B° — > K®ir° decay no charged particles originate from 
the decay vertex, we compute the decay point of the B rcc 
using the trajectory, obtained from the reconstructed 
ATg decay vertex and momentum vector, and the average 
e + e~ interaction point [38]. We have verified that all 
correlations between the fit variables are negligible, and 
so construct the likelihood function as a product of one- 
dimensional PDFs. Residual correlations are taken into 
account in the systematic uncertainty, as explained be- 
low. 

The PDFs for signal events are parameterized based 
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on a large sample of fully-reconstructed B decays in data 
and from simulated events. For background PDFs, we 



take the functional form from the background-dominated 
sideband regions in the data. The likelihood function is: 



J 



£(S K o o,C K o o,N s ,N B ,fs )/!'«) = 



,-(N s +N B ) 

AM 



(9) 



i£ 9 



r 



where the N selected events are partitioned into two 
subsets: the index i £ g indicates events that have Ai 
information, while i € b events do not have At infor- 
mation. Here, /| (/g) is the fraction of signal (back- 
ground) events that are in the subset g, and fs = l—fi 
(/ B = 1 — /b) are the corresponding signal (background) 
fractions in the subset b. The parameter Ns (Nb) is the 
number of signal (background) events. The probabili- 
ties Ps and Vb are products of PDFs for the signal and 
background hypotheses evaluated for the measurements 
Xi = {tub, L 2 /L o, cos9 B , flavor tag, tagging 

category} and yi — {At, a&t}- The corresponding PDFs 
for events without At information are V'§ and P B . De- 
tailed descriptions of Vs- V B , V' s . and V' B are given in 
Ref. 1. The vector a represents the set of parameters 
that define the shapes of the PDFs. Along with the CP 
asymmetries S K o^o and C K o n o, the fit extracts the yields 
Ns and N B , the fraction of events /| and f B , and the pa- 
rameters of the background PDFs. 



IV. RESULTS AND SYSTEMATIC 
UNCERTAINTIES 

. B° TT+TT" and B° 



TABLE II: Results for the B° -> h + h'~ decay modes. Un- 
certainties on the signal yields N s i s are statistical. For the 
CP-violation parameters, the first uncertainties are statisti- 
cal, and the second are systematic. 



Mode 




CP-violation parameters 




7T + 7T 


1394 ± 54 


S v+V - = -0.68 ± 0.10 ± 0.03 
C n+7I - = -0.25 ± 0.08 ± 0.02 


B l> — > 


K + 1T~ 


5410 ± 90 


A K - V+ =-0.107 ±0.016™ 




K + K- 


7 ±17 





The event yields and CP-violation parameters are 
listed in Table [TTJ The correlation coefficient between 



and CL 



is found to be 



relation between 0^+^- and Ak~ 



-0.056, and the cor- 
- is 0.019. We show 



B — » Kir, and qq background in Fig. [2j where the 
s Plots [39J weighting and background-subtraction tech- 
nique is used to display a distribution for a particular 
type of event. The direct CP asymmetry in B° — > K + w~ 
is apparent in the AE distributions, plotted separately 
for B° and B° decays in Fig. [3J We show the distri- 
butions of At for B° — > K^n 11 signal and background 
decays in Fig. [4] In Fig. [5j we show the distribution of 
At separately for B — » ii + it~ events tagged as B° or 
B°, as well as the asymmetry a(At) of Eq. |T]). The re- 
sults for SV+tt- and C v + V - are shown in Fig. |6j along 
with confidence-level contours corresponding to statisti- 
cal significances ranging from la to 7a. Our measure- 
ment excludes the absence of CP violation in B° — » ir + n~ 
(S n + n - = 0, C^+k- = 0) at a confidence level corre- 
sponding to 6.7a, including systematic uncertainties. 



TABLE III: Summary of systematic uncertainties on A K - 7T + . 
To address the Ak-^+ bias due to hadronic interactions of 
charged kaons with the detector material, we shift the A K -^+ 
value obtained in the fit by +0.005. 



Source 


A K ~7T + 


Material interactions 


+0.005 -0.003 


B c and dE/dx PDFs 


0.002 


Alternative DIRC parameterization 


0.002 


Potential bias 


0.001 


Total 


+0.006 -0.004 



Systematic uncertainties for the direct CP asymmetry 
Ak--k+ are listed in Table [TTT1 Here, Ak--k+ is the fit- 
ted value of the event-yield asymmetry _4 I , aw 



+0.005 
-0.003 



to account for a bias that arises 



shifted by +0.005 
from the difference between the cross sections of 
and K~ hadronic interactions within the BABAR detec- 
tor. We determine this bias from the MC. The bias is 
independently verified with a calculation based on the 
known material composition of the BaBAR detector [23| 
and the cross sections and material properties tabulated 
in Ref. (29|. The corrected K^-k^- event-yield asymme- 
try in the background, where no observable CP violation 
is expected, is -0.005 ± 0.004 (stat)!";^ (syst), consis- 
tent with zero. Uncertainties on the 9c and dE/dx dis- 



the rriES, AE, and T distribution for the B — > 7T7t, tributions are obtained from the D° -+ K 7r + control 
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FIG. 2: sVlots of the (left column) rriES, (center column) AE, and (right column) Fisher discriminant T distributions for 
(top row) B° —¥ it + tt~, (middle row) B° — > K + ir~ , and (bottom row) qq background candidates. The points with error bars 
show the data, and the lines represent the PDFs used in the fit and reflect the fit result. The structure to the left of the 
signal AE peak for B° — > tt + it~ is consistent with the expected background from other charmless modes, which is negligible 
for AE > —0.10 GeV. In the calculation of AE for B° — > K + ir~ , the kaon candidate is assigned the pion mass. 




-0.1 0.1 

AE (GeV) 

FIG. 3: s Plots of the AE distribution for signal K^-k^ 
events, comparing (blue solid lines, filled circles) B° and (red 
dashed lines, empty circles) B° decays. The points with error 
bars show the data, and the lines represent the PDFs used in 
the fits and reflect the results of the fits. 



TABLE IV: Summary of systematic uncertainties on S'^+ 7r - 
and C^ +7r -. 



Source 






DIRC 9c 


0.0064 


0.0050 


DCH dE/dx 


0.0032 


0.0037 


Signal At 


0.0199 


0.0055 


SVT local alignment 


0.0004 


0.0002 


Boost/detector z size 


0.0021 


0.0013 


PEP-II beam spot 


0.0028 


0.0014 


B flavor tagging 


0.0146 


0.0138 


Am d , t b o [29] 


0.0004 


0.0017 


Potential bias 


0.0041 


0.0043 


CP violation in Stag decays 


0.007 


0.016 


Total 


0.027 


0.023 



sample, and contribute 0.002 to the systematic uncer- 
tainty on Ak-tt+ ■ An additional uncertainty of the same 
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FIG. 4: s Vlots of the At distribution for (top) signal K ± ir^ 
and (bottom) background events. The points with error bars 
show the data, and the lines represent the PDFs used in the 
fit and reflect the fit result. 



magnitude is obtained by adding a bifgurcated-Gaussian 
component to the two-Gaussian 6q PDF. We use a com- 
bination of MC events and parameterized experiments to 
test for a potential bias in the fit, for which we estimate 
an uncertainty of 0.001. 




FIG. 5: aVlots of the At distributions for signal tt + it~ events 
tagged as (top) B° or (middle) B°, and (bottom) their asym- 
metry a(At), from Eq. fT}. The points with error bars show 
the data, and the lines represent the PDFs used in the fit and 
reflect the fit result. 



Systematic uncertainties for the CP asymmetries 
SjY+Tf— and C~.+~.— are listed in Table IIVI The largest 
uncertainties on S^+^- are due to the At and £?-flavor- 
tagging parameters, and are determined by varying the 
At resolution function parameters and the flavor-tagging 
parameters by their uncertainties. The largest C„+„- un- 
certainty is due to the effect of CP violation in the Stag 
decays 40]. The effect of SVT misalignment is deter- 
mined by reconstructing events with shifted alignment 
parameters, and the uncertainties due to the machine 
boost and detector size are obtained by scaling At by 
1.0046. We evaluate uncertainties due to the measure- 
ment of the beam spot by shifting its position in the 
vertical direction by 20 //to, and those due to the knowl- 
edge of the B° — B° mixing frequency and the B° lifetime 
are determined by varying these parameters within their 
uncertainties [2fj| . The uncertainties due to particle iden- 
tification and potential fit bias are evaluated as described 
above for Ak-k+- 



B. B° ttV Results 

Results from the ML fit for the B° — > tt tt decay 
mode are summarized in Table fVl s Vlots of toes, AE, 
and NN for B° — > 7r°7r° are shown in Fig. [7J and for the 
qq background in Fig. [5] 

The various systematic uncertainties for the B° — > 
7r°7r° decay mode are listed in Tables I VI I and I VIII The 
uncertainty in the efficiency is dominated by a 3% sys- 
tematic uncertainty per 7r°, estimated from a study of 
t — > tvk Vt decays. An uncertainty of 1.0% is due to 
the resolution of the signal shape, and an additional un- 
certainty of 0.5% is due to the limited knowledge of the 
toes and AE peak positions in data. These are esti- 
mated by shifting the toes an d AE means and resolu- 
tions by amounts determined from MC-data comparison 
in a control sample of B + — > 7r + 7r° events. An uncer- 
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TABLE V: Results for the B — > 7r°7r and B — > K"tv decay modes, showing the signal yield iV s ig, efficiency, branching 
fraction, and CP-violation parameter C for each mode. When two uncertainties are given, the first is statistical and the second 
is systematic. Uncertainties for the signal yields are statistical, and those for the efficiencies are systematic. 







Efficiency (%) Branching fraction (10 b ) 


C 


B l> - 


-> ttV" 247 ± 29 


28.8 ±1.8 1.83 ±0.21 ±0.13 


-0.43 ±0.26 ±0.05 


B°- 


* K° s -k° 556 ± 32 


34.2 ±1.2 5.1 ±0.3 ±0.2 
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FIG. 6: S^+^- and C w + n - in B — > tv + tt~ decays, show- 
ing the central values (point with error bars) and statistical 
confidence-level (C.L.) contours for 1 — C.L. = 0.317 (la), 
4.55xl0" 2 (2a), 2.70xl0" 3 (3a), 6.33xl0" 5 (4a), 5.73xl0~ 7 
(5a), 1.97 x 10~~ 9 (6a) and 2.56 x 10~ 12 (7a), calculated from 
the square root of the change in the value of — 21n£ with re- 
spect to its value at the minimum. The unit circle represents 
the physical region S 2 +7r _ ± C 2 +^_ < 1. 



tainty of 1.5%, determined from the -Bfl av sample, is due 
to the |cos# s | requirement. A 1.1% uncertainty is as- 
signed to the number of BB events in the data sample. 
Systematic uncertainties involving the ML fit are eval- 
uated by varying the PDF parameters and refitting the 
data. These contribute an uncertainty of 8.3 events to 
the branching-fraction measurement and an uncertainty 
of 0.055 to CLo_o. 




o 2 4 6 8 i( 

NNbin 



FIG. 7: B — > it 7T signal plots with background subtracted 
using the s Vlots technique. From top to bottom: rriES, AE, 
and AW. The points with error bars show the data, and the 
line in each plot shows the corresponding PDF. 



C. B° K%-ir° Results 

The efficiency and branching fraction measured for the 
B° — » K®tt° decay mode are summarized in Tablc|V] (CP- 
violation parameters have been reported in Ref. j4|). 

We show s Vlots of m m i ss , m^, L2/L0, and cosd* B 
for signal events in Fig. |H] and for background events in 
Fig. [101 

The systematic uncertainties on the branching frac- 



tion B(B° — > K®ir°) are summarized in Table IVIIII 
The uncertainty on the efficiency of the K s reconstruc- 
tion is obtained from detailed comparison of inclusive 
Kg candidates in data and MC. The ir° efficiency un- 
certainty is evaluated from the ratio of branching frac- 
tions B{D° K-Tr+TT )/B{D° -> K~-k + ). To compute 
the systematic uncertainty associated with the statistical 
precision on the parameters of the likelihood function, we 
shift each parameter by its associated uncertainty and 
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background plots with signal subtracted 



FIG. 8: 5° 

using the s Vlots technique. From top to bottom: m^s, AE, 
and NN. The points with error bars show the data, and the 
line in each plot shows the corresponding PDF. 



7r°7r° sig- 



TABLE VI: Systematic uncertainties on the B° 
nal yield N v o n o and direct CP asymmetry C^o^o. The total 
uncertainty is the sum in quadrature of the individual uncer- 
tainties. 



Source 






Peaking background 


4.9 


0.030 


Tagging 


0.35 


0.034 


Background shape 


5.5 


0.023 


Signal shape 


3.8 


0.020 


Total fit systematic uncertainty 


8.3 


0.055 



TABLE VII: Relative systematic uncertainties on the B° — > 
branching fraction. The total uncertainty is the sum 
in quadrature of the relative uncertainties on the signal yield 
(from Table |VT)) . the signal efficiency, and the number of BB 
pairs. 



Source 


B(B U ->■ VY J ) 


Signal yield syst. uncertainty 


3.4% 


7r u efficiency 


6.0% 


| cos 8s | selection 


1.5% 


neutrals resolution 


1.0% 


jd.es and AE shape 


0.5% 


Number of BB pairs 


1.1% 


Total systematic uncertainty 


7.2% 




Sioo 
w 



5.15 5.2 5.25 5.3 
m miss (GeV/c 2 ) 



5.2 5.3 5.4 
m B (GeV/c 2 ) 




repeat the fit. For At and the tagging parameters, the 
uncertainty is obtained from the fit to the -Bfw sample, 
while for the other parameters it is obtained from MC. 
This uncertainty accounts for the size of the sample used 
for determining the shape of the likelihood function in 
Eq. @. A systematic uncertainty associated with the 
data-MC agreement in the shape of the signal PDFs is 



FIG. 9: s Plots of the (a) m m i ss , (b) tub, (c) L2/L0, 
(d) cos #3, and (e) At distributions for signal events in the 
B° — > KsTv sample. The points with error bars represent 
the data, and the lines show the shapes of signal PDFs as 
obtained from the ML fit. 
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cay channels that dominate this type of background, no- 
tably B+ -> p + K° s , B+ -> K*+n°, and B+ -> Ar°7r°7r+. 
We embed these simulated B-background events in the 
data set and find the average shift in the fit signal yield 
to be +5.2 events. We adjust the signal yield accordingly 
and use half of the bias as a systematic uncertainty. 

For the branching fraction, additional systematic un- 
certainties originate from the uncertainty on the selection 
efficiency, the number of BB pairs in the data sample 
(1.1%), and the branching fractions B(K° S — » 7r + 7r~) and 
B(tt ->77) H3. 



TABLE VIII: Summary of dominant contributions to the sys- 
tematic uncertainty on the measurement of B(B° — > Kgir ) 



Source 


a(B(B" A> u )) (%) 


7r ,J efficiency 


3.0 


K® efficiency 


0.5 


Selection criteria 


1.5 


PDF-parameters precision 


0.22 


Shape of signal PDFs 


0.45 


BB background 


0.47 


Correiations 


0.40 


Resolution function 


0.49 


Number of BB pairs 


1.1 


Total 


3.7 



FIG. 10: s Tlots of the ( (b) m B , (c) L 2 /L , 

(d) cos 9%, and (e) At distributions for background events 
in the B° — > K®-k° sample. The points with error bars repre- 
sent the data, and the lines show the shapes of signal PDFs 
as obtained from the ML fit. 



evaluated by taking the largest deviation observed when 
the parameters of the individual signal PDFs for m m i ss , 
ttib, L2/L0, and cos0|j are allowed to vary in the fit. 
The output values of the PDF parameters are also used 
to assign a systematic uncertainty to the efficiency of 
the event selection requirements on the likelihood vari- 
ables, by comparing the efficiency in data to that in the 
MC. We evaluate the systematic uncertainty due to the 
neglected correlations among fit variables using a set of 
MC experiments, in which we embed signal events from 
a full detector simulation with events generated from the 
background PDFs. Since the shifts are small and only 
marginally significant, we use the average relative shift 
in the yield as the associated systematic uncertainty. 

In the fit we neglect background from B decays, which 
is estimated from simulation to contribute of order 0.1% 
of the total background. To account for a bias due to this, 
we study in detail the effect of a number of specific B de- 



V. RESULTS FOR Aa^ AND a 

We combine our results for B(B° — > tt°tt ) with the 
branching fractions B(B° -> tt+tt") = (5.5 ± 0.4 ± 0.3) x 
10" 6 and B{B ± ->■ 7 r ± 7r°) = (5.02 ± 0.46 ± 0.29) x 10~ 6 
previously measured by BABAR 0, [l5[ to evaluate the con- 
straints on both the penguin contribution to a and on the 
CKM angle a itself. Constraints are evaluated by scan- 
ning the parameters lAo^l and a, and then calculating 
the x 2 for the five amplitudes (A +a , A + ~ , A 00 , A+~ , 
A 00 ) from our measurements and the isospin-triangle re- 
lations 10]. Each x 2 value is converted to a confidence 
level, shown in Fig. [TT] for Aa M and a. The a plot ex- 
hibits six clear peaks, a result of the eight-fold trigono- 
metric ambiguity in the extraction of a and the fact that 
two pairs of peaks are nearly merged. The upper bound 
on IAq^I is 43° at the 90% C.L., and the range [23°, 67°] 
in a is excluded at the 90% C.L. The point a = 0, which 
corresponds to no CP violation, and the values of a near 
or 7T can be excluded with additional physics input I6l 4l|. 
If we consider only the solution preferred in the SM [42| , 
a lies in the range [71°, 109°] at the 68% C.L. This is 
consistent with the more restrictive constraints on a ob- 
tained from analysis of the B — > pp system I43| . as well 
as those from B° -t (p7r)° [il] and B° -> a-nr [45l]. 



17 




'0 50 100 150 
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FIG. 11: (Top) Constraint on Aa^^ = a — a e B> expressed as 
one minus the confidence level as a function of |Aa W7r |. We 
find an upper bound on |Aa ww | of 43° at the 90% C.L. (Bot- 
tom) constraint on the CKM angle a. We exclude the range 
[23°, 67°] in a at the 90% C.L. Only the isospin-triangle rela- 
tions and the expressions in Eq. ([TJ are used in this constraint. 

VI. CONCLUSIONS 

We measure the CP-asymmetry parameters 

S^- = -0.68 ±0.10 ±0.03, 
CV +7r - = -0.25 ±0.08 ±0.02, 
Ak-«+ =-0.107 ±0.016™, 
CWo = -0.43 ±0.26 ±0.05 

and CP-averaged branching fractions 

B(B° -> 7r°7r°) = (1.83 ± 0.21 ± 0.13) x 10~ 6 , 
B(B° -> A'°7r°) = (10.1 ± 0.6 ± 0.4) x 10~ 6 . 

We find a 68% C.L. region for a of [71°, 109°] and ex- 
clude values in the range [23°, 67°] at the 90% C.L. We 



observe direct CP violation in B° — > K + tt~ with a signif- 
icance of 6.1cr and in B° — >• 7r + 7r~ with a significance of 
6.7cr, including systematic uncertainties. Ignoring color- 
suppressed tree amplitudes, the charge asymmetries in 
K + n~ and K + tt° should be equal [21) . which is not 
supported by recent BABAR and Belle data 0, H H^. 
These results might indicate a large color-suppressed am- 
plitude, an enhanced electroweak penguin, or possibly 
new-physics effects [471 ] . 

Our result for B(B° — > K°tt°) is consistent 
with the sum-rule prediction [U |22( i3(i^r°7r°) sr = 

§ (B(K+n-) + a. [B(K°7T+) - 2B(K+tt )}]^ = (8.4 ± 
0.8) x 10" 6 , obtained using the currently published re- 
sults [1, dB-Ell for the three B —> Kn rates on the right- 
hand side of this equation and the lifetimes r + and To of 
the charged and neutral B mesons. 

The resultspresented here supersede those of our prior 
publications 
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